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For the Gaussian and Laguerre random matrix ensembles, the probabihty density 
function (p.d.f.) for the Unear statistic Y^f=i {xj — {x)^ is computed exactly and 
shown to satisfy a central limit theorem as — > cxo. For the circular random matrix 
ensemble the p.d.f.'s for the linear statistics ^ J2jLii^j~'^) ^^^^ ~ J2jLi log 2| sin 6*^/21 
are calculated exactly by using a constant term identity from the theory of the Selberg 
integral, and are also shown to satisfy a central limit theorem as — > oo. 
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1 INTRODUCTION 

A number of rigorous results have recently been established regarding probability density func- 
tions and associated fluctuation formulas for linear statistics in random matrix ensembles. We 
recall that with Ai, . . . , Atv denoting the eigenvalues of a random matrix, a linear statistic is any 
stochastic function A which can be written in the form A = ^jLi a{Xj). Corresponding to the 
statistic A is the probability function (p.d.f) 



N .oo poo f ^ \ 

P{u):=l[ dXi... dXNS\u-Y,a{Xj)]w{Xi,...,XN), 



fi.i) 



where W denotes the p.d.f. for the eigenvalue distribution, for the event that u = J2jLi 

An evaluation of the variance of P{u) is referred to as a fluctuation formula, although this term 

is sometimes used to refer to an evaluation of P{u) itself. 

When a{x) = X[-l/2,l/2] - {X[-l/2,l/2]) , where ( ) denotes the mean and xb = 1 if a: G 5, 
Xb = otherwise, we have that the linear statistic A represents the deviation in the number 
of eigenvalues from the mean number in the interval [—1/2,1/2]. For N x N Gaussian random 
matrices, scaled so that the mean eigenvalue spacing in the bulk of the spectrum is unity, it has 
been proved by Costin and Lebowitz [1] that 

P ((^logl) ' ~ e-"'/2 as / ^ cx), (1.2) 



where /? = 1,2 and 4 for random symmetric, Hermitian and self-dual real quaternion matrices 
respectively. Before scaling, the eigenvalue p.d.f. for Gaussian random matrices is proportional 
to 

N 

n^-''''/' n iA.-A,r. (1.3) 

l=\ l<j<k<N 
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This can be interpreted as the Boltzmann factor for a log-gas system in equihbrium at inverse 
temperature (3. Prom this interpretation, heuristic arguments based on macroscopic electrostat- 
ics have been devised [2,3] which support the validity of (1.2) for all /? > 0. 

Another, more general class of rigorous results for the evaluation of (1.1) for random matrix 
ensembles has been obtained by Johansson [4,5] (see also [6] in the case (5 = 2). In these results, 
instead of scaling the eigenvalues so that the mean spacing is unity, the scale is chosen so that 
ss N ^ CO the support of the density of the eigenvalues is the finite interval (—1, 1). Since for 
large the density of eigenvalues p{x) implied by (1.3) is given by the Wigner semi-circle law 
(see e.g. ref. [7]) 

this is achieved by the scaling Xj i— > (2iV)^^^Aj in (1.3). More generally, with the exponent in 
(1.3) replaced by an even-degree polynomial which is positive for large A;, a scale can always be 
chosen so that the support is the finite interval (—1, 1) [4]. Moreover, independent of the details 
of the polynomial, Johansson proved that the p.d.f. (1.1) in the scaled N ^ oo limit tends to a 
Gaussian: 

support p(a:) = ( — 1,1) 

where 



lim Piu) = e-" /^'^ (1.5a) 

AT— >oo 
support p{x) = { — l,l) 



^ =1^ I ^-^'t: / ^ z—. (i-5b) 



Pn^ J-i ' (1 -x2)i/2 x-y 

provided cj^ is finite. The value of the variance was predicted in earlier work due to Brezin and 
Zee [8] in the case /3 = 2, and Beenakker [9] for general /3 > 0. 

Both (1.2) and (1.5) can be interpreted as central limit theorems. However, unlike the 
classical result the standard deviation is no longer proportional to ^/N (or \/I in the case of 
(1.2)). Rather the fiuctuations are strongly suppressed, being independent of A'" altogether in 
the case of (1.5). 

It is the purpose of this work to extend the rigorous results relating to the p.d.f. (1.1) for 

random matrix ensembles. Attention will be focusscd on two types of linear statistics: a{x) = x 
and a{x) = — log \x\, which both arise naturally within the log-gas interpretation of (1.3). Indeed 
a(x) = X corresponds to the dipole moment while a{x) = — log |x| corresponds to the potential 
at the origin. This latter statistic has been considered for the two-dimensional one-component 
plasma (i.e. the two-dimensional generalization of the log-gas) by Alastuey and Jancovici [10]. 

In Section 2 we will show how (1.1) with W proportional to (1.3) can be computed exactly 
for a{x) = X. Also, we calculate (1-1) exactly for a{x) = x — [x) with W proportional to 

N 

J|^/5a/2^-^../2 -Q \x,-x,f, (1.6) 
1=1 l<j<k<N 

Xj > {j = 1, . . . , N), corresponding to the Laguerre random matrix ensemble. In the scaled 
AT — > oc limit, the result (1.5) with a{x) = x is reclaimed in both cases. 
In Section 3 we consider (1.1) with W proportional to 

n {2\sm{e,-ej)/2\f, (1.7) 

l<j<k<N 

< Oj < 27r (j = 1,...,A^), corresponding to Dyson's circular ensemble of unitary random 

matrices. In the case (3 = 2 the large- A^ behaviour of P{u) for a{9) = {0 — tt)/2 and a{9) = 
— log 2 1 sin 0/2 1 is computed using the so-called Fisher-Hartwig conjecture [11] (which is now a 



2 



theorem [12,13]) from the theory of Toeplitz determinants. For general f3, P{u) is evaluated 
for these statistics by using a constant term identity due to Morris [14] which is equivalent to 
the well known Selberg integral [15]. We conclude in Section 4 with an interpretation of the 
variance of the dipole moment statistic calculated in Section 2 as a susceptibility in macroscopic 
electrostatics. 

2 DIPOLE MOMENT STATISTIC IN THE GAUSSIAN AND 
LAGUERRE ENSEMBLE 

Rather than study the p.d.f. (1.1) directly, we consider instead its Fourier transform 

/oo 
-OO 

/OO roo 
dxi e^'^^^^i) . . . / dxN e**^"(^^) W^(xi, . . . , xn). (2.1) 
-oo J —oo 

Note that (2.1) can be interpreted as the canonical average of the Boltzmann factor for a one- 
body external potential ika{x)/f5. 

2.1 Gaussian random matrices 

There are three distinct random matrix ensembles, in which the matrices X are real symmetric 
{13 = 1), Hermitian {(3 = 2) and self-dual quaternion real {(3 = 4), and the joint distribution for 
their elements is proportional to e~^'^^'^'^^/'^ . The corresponding eigenvalue p.d.f. is proportional 
to (1.3). Scaling the eigenvalues \j ^ V 2NXj so that the support of the density is (—1, 1), we 
have from (1.3) and (2.1) that for the linear statistic corresponding to a{x) = x 



^ ^-"^ l<j<k<N 

where the normalization C is such that P{0) = 1. 

Now for /3 even the integrand in (2.2) is analytic and decays sufficiently fast at infinity that 
we can shift the contours of integration from the real line to Xj = ik/2f3N + tj (j = 1, . . . ,N), 
which gives that the integral is independent of k. Since the integral in (2.2) divided by C is a 
bounded analytic function of /? for Re(/3) > 0, it follows from Carlson's theorem [16] that the 
integral is independent of k for all Re(/3) > 0. Thus for all /3 > and each = 1, 2, ... we have 
the simple result that 

P{k) = e-^^l^^ (2.3a) 

and so 

= [^^"^ e-P^\ (2.3b) 

Since (2.3b) is independent of N it trivially remains valid in the N ^ oo limit, and agrees with 
(1.5) provided = 1/2/3. 

The value of cr^ is easily computed from (1.5b) by recalling (see e.g. [17]) that the solution 
of the integral equation 

/■I Mil) 

(2.4) 



x-y 
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is = — y'^. Thus with a{x) = x, substituting (2.4) in (1.5b) we have 
as required. 

In Fig. 1 we ihustrate (2.3b) by empiricaUy calculating P{u) for the eigenvalues of 5,000 2x2 
matrices from the Gaussian Orthogonal Ensemble of random real symmetric matrices (/3 = 1) 
and comparing the empirical p.d.f. to the theoretical value, eq. (2.3b) with (3 = 1. 

2.2 Laguerre random matrix ensemble 

If X is a random nx m (n > m) matrix with Gaussian entries, which are either all real (/? = 1), 
complex (/3 = 2) or real quaternion (/3 = 4) with joint distribution of the elements proportional 
to e~''^'^(^^^)/^, then the eigenvalue p.d.f. of X^X is proportional to (1.3) with N = m and 
a = n — m + x/3 {xp = ~1) Oj 1 for /3 = 1, 2, 4 respectively). For this so called Laguerre ensemble 
it is known (see e.g. [6]) that the support of the eigenvalue density is (0,4A^). Scaling the 
eigenvalues Xj i-^ 2NXj, so that the support is (0,2) (the important point here is that the length 
of the interval is 2, as is required for the validity of (1.5b)), (2.1) with a{x) = x — {x) reads 



1 /* OO P oo 

P{k) = ^ dxixf /'e-^^-i+^^(-i-<-» . . . / cZxjv4"/'e-^^--+^'=(---<-» ff \xk - x,f , 

(2.6) 

where C is such that P(0) = 1 and 
(x) := i r dx, /^e-/^^- ...r dxM x%^'\-^^^- f 1 x\ J] 1^^ " ^if ■ (2-7) 
To evaluate (2.7), we note that a simple change of variables gives 

poo poo 

/ dxi e-'^^^i . . . / dxNX%e-i'^''^ TT |xfe - Xjf 
•^0 -^0 i<j<ik<Ar 

rOO rOO 

^ ^-N{a+l)-m(N-l)/2 dxiX^-'''''... dXivX^e-^^^ TT |Xfc - X,- 1^ (2.8) 

^0 Jo 



l<j<k<N 



and that the integral in (2.7) can be obtained from the left hand side of (2.8) by differentiation 
with respect to fi and setting a = Pa/2, b = fiN and fi = 1. Applying the same operation to 
the right hand side of (2.8) gives 

(x) = i + -i- f a - 1 + I V (2.9) 



2 2N \ P 

Now, the non-trivial dependence on k in (2.6) occurs only in the factors e~('^^+*^)^J . This 
dependence can be taken outside the integral by the change of variables {(3N — ik)xj = PNyj 
(this operation is immediately valid for P even; it remains valid for Re(/3) > by Garlson's 
theorem) to give 

Taking the inverse transform we have 

P{u) = r(^i^2^^)) (W''^^H^ + iV(x))^^^<^)-^e-/^^("+^<-)), (2.11) 
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for u > —N{x), P{u) = for M < ■~N{x). In the limit ^ oo we see that (2.10) tends to 
(2.3a) as expected, thus explicitly demonstrating the universality feature of (1.5). 

In Fig. 2 we illustrate (2.11) by numerically calculating a histogram for the p.d.f. oiJ2jLii^j~ 
(A)) for 3,000 matrices X'^X, where X is a 3 X 2 real rcctanglar matrix with entries chosen with 
p.d.f. -^e"^^*^^, and comparing it against the theoretical prediction (2.11) with P = 1, N = 2 
and (x) = 1/2 + 1/2N. 



3 DIPOLE MOMENT AND POTENTIAL STATISTIC IN THE 
CIRCULAR ENSEMBLE 

3.1 The Vciriance 

For a general linear statistic A, 

_2 . /r A /4\\2\ 



{(A - {A)Y 

ji.jiyai.)a(y){f,H.,y)+,,(m^-y)). (3.1) 



where pi{x) denotes the density, the truncated two particle distribution and / the allowed 
domain for the particles (eigenvalues). For the circular ensemble we have I = [0, 27r) and due to 
the periodicity we can write 

oo 

pl{x,y)+pi{x)5{x-y)= rj^'"'^>. (3.2) 

n=— oo 

Substituting (3.2) in (3.1) gives 

oo oo 

(72 = 27r ^ |a„|V„, where a(x) = ^ a„e'"^, a;G[0,27r). (3.3) 

n=— oo n=— oo 

Consider now the particular linear statistics with a\{x) = ^{x—tt) and a2{x) = — log 2| sina;/2|. 
We have 

sin ma; cosmx 
ai{x) = 2^ —— and a2{x) = }^ , (3.4) 

m=l m=l 

SO to analyze cr^ from the formula in (3.3) it remains to specify the behavior of r„. This can be 
done using an heuristic electrostatic argument combined with linear response theory (see e.g. 
[18]), which gives ^ I'^I/tt/J for < |?7.| ^ 0(iV). Substituting in (3.4) we therefore have 

cN ^ . 

for both the dipole moment and potential statistics. Linear response arguments [19,20] also give 
the prediction that in an appropriate macroscopic limit (which here corresponds to ^ oo) 
the p.d.f. of any linear statistic will be Gaussian, thus suggesting that for the dipole moment 
and potential statistics 

P (-logiV) u ~( — ) e-"'/2 as N^oo. (3.6) 



p J ; \2ttJ 
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3.2 The case p = 2 

^From the Vandermonde determinant expansion 



det 



^t%,-, ^= n (-^--.) (3-7) 

it is straightforward to show (see e.g. [21]) that with W given by (1.7) and (3 = 2 (2.1) can be 
rewritten in terms of a Toephtz determinant: 

(•27r 



P{k) = Dm 





:= det 


\-[ 






2ttJo 



j,k=l,...,N 



(3.8) 



Using an asymptotic formula for the large- form of D^, known as the Fisher-Hartwig conjec- 
ture [11-13], we can rigorously establish (3.6). 

Proposition 3.1 (Fisher-Hartwig conjecture) In (3.8) let 

R R 
ika{9) = g{e) - i^hriir - {9 - Or)) + ^ar log |2 - 2 cos(^ - ^^)| (3.9) 



r=l 



r=l 
|2 



and assume g{9) = J2'^-oo9p^ where Z^^-ooblbpl < oo- Then for Re(Q!r.) > —1/2 and 



D 



N 



Jka{9) 



:5oA'eEf=i(«'-''')i°g^_E 



(3.10) 



where E is independent of N. To specify E, write g{9) — go = g'+(e*^) -f- g-{e *^) where g+{z) 
E^i 9pZ^ and g-{z) = Epi-oo 9pZ^- Then 

R 

E = e^kLi^9k9-k Jl g-(ar.+6r)g-(e'^'-)g-(ar-6r)3+(e-'^) 



r=l 



|-(ar+&r)(as-&s) G{1 + Gr + br)G{l + — br) 
l<r<s<R r=l 

where G is the Barnes G-function defined by 



G{z + 1) = (27r)^/' exp (-z/2 - (7 + l)zV2) H ( ^ + Z ) exp(-z + z^/k) 



k=l 



(3.12) 



(7 denotes Euler's constant), which has the special values G{1) = G{2) = 1 and satisfies the 
functional relation 

G{z + 1) = r{z)G{z). (3.13) 

First consider the application of (3.9) to the calculation of the p.d.f. for the potential statistic. 
From (3.8) we have 

(3.14) 



P{k)=DN [(2 -2 cos 0)^^=/^' 



so we take g{e) = 0, i? = 1, 61 = 0, ai = f (^logAT) ' and 0i = in (3.9) and (3.10) to 
conclude that 

k 



(ilogiV) 



1/2 



~ e 



(3.15) 



which is equivalent to (3.6) with f3 = 2. 

For the dipole moment statistic we have from (3.8) that 



P{k) = Dn 



(3.16) 



Thus we take g{e) = 0, i? = 1, ai = 0, 61 = | and 6'i = in (3.9) and (3.10) to conclude that 
(3.15) again holds, as predicted from (3.6). 
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3.3 General (3 

Here we will show that P{k) for the dipole moment and potential statistics in the circular 
ensemble can be given in closed form for general /?, and prove that the expected asymptotic 
behaviur (3.6) is valid for general rational (3 (at least). Our chief tool for this purpose is a 
constant term identity of Morris [14] , written as the multidimensional integral evaluation [22] 

M„(a, 6, c) := n — deie^^'^-'^'^/^\l + e^'^ f^' [] 1^*'' " e^'^' l'' 

r(c(j + l) + l)r(rj+a + 6+l) 
/J, V{cj + a + l)V{cj + h + l)r(l + c) ■ ^ - > 



First consider the dipole moment statistic. Prom (1.7) and (2.1) we have 



Z=l'^° l<j<k<N 

MN{k/2,-k/2,f3/2) (r(/3i/2 + 1))^ 

M^r(o, 0, p/2) r(/3j72 + fc/2 + i)r{pj/2 -k/2 + iy ^ ■ > 

Although (3.18) is a closed form evaluation, it is not convenient for the determination of the 
large-A^ asymptotics. This same problem has been faced in an earlier application of the Morris 
constant term identity [23] . Its resolution is to make use of the identity 

fr' r(a + 1 + jc) = c'^NiN-l)/2+N(a+l/2) (^2^yNic-l)/2 TT GiN+ {a + l)/c) 

valid for c a positive integer. We can use (3.19) in (3.18) for all rational /3, (3/2 = s/r with s 
and r relatively prime positive integers say. This is done by replacing A'' by rN and noting 

rN-l , ■ . r-17V-l , . 

n r(«+i+f) = nnr(«+i+f (3.20) 

j=0 ^ ^ p=0 fe=0 ^ ^ 

Thus for /3/2 = s/r and with N replaced by riV, (3.18) can be rewritten to read 

r(/:) TT fr ^^^^ ^ ^ + + + - fe/2)/g) .3 2l^ 

/ioti G"!^ + + + fc/2)/s)G'(A +^)s/r + (^ - A;/2)/s)(G(l/s +j)s/r))2' 

The large-TV asymptotics of -P(A;) can be deduced from (3.21) by making use of the asymptotic 
formula [24] 

( g(^ + 6+!) ) ~ (& - + ^ log 2^ + ((« - + log + o(l). (3.22) 

This gives 

log . 4 log^ + V V log gWr + (i + «./2)A)GWr + (i-V2)/.) 



which implies 

1/2 



(^logiv) ; 
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in precise agreement with the Fourier transform of the expected asymptotic behaviour (3.6). 

Let us now turn our attention to the potential statistic. In this case, from (2.1), (1.7) and 
(3.17) we have 

1=1 "'^ l<j<k<N 

MNi-ik/2, -ik/2, (3/2) ^ T{f3j/2 - ik + l)T{f3j/2 + 1) 

M^(0,0,/3/2) jj^ {r{Pj/2-ik/2 + l)f ^' ' 

Proceeding as above, for rational P {P/2 = s/r), by using (3.19) and (3.20) we can rewrite (3.25) 
as 

_ TT fr G{N + ps/r ik)/s)G{N + ps/r + l/s){G{ps/r + {I - ik/2)/s)f 

piozii {G{N +ps/r + {l-ik/2)/s)fG{ps/r + {l-ik)/s)G{ps/r + l/s) 

The asymptotic formula (3.22) now gives that for A'" — oo 

P^oi^i G(jps/r + {l-ik)/s)G(jps/r + 1/8) 

which again implies the anticipated result (3.24). 

In Fig. 3 we illustrate (3.18) by numerically constructing a histogram for the p.d.f. of 

\Yl^=i{Oj — vr) for 5,000 3 x 3 matrices from the CUE (these are constructed according to 
the procedure specified in [25]), and comparing it against the inverse transform of (3.18) in the 
case N = 2,, 13 = 2. 

4 RELATIONSHIP BETWEEN THE VARIANCE OF THE 
DIPOLE MOMENT STATISTIC AND THE SUSCEPTIBIL- 
ITY 

In this final section we will show how the result (2.5) can be anticipated from its interpretation 
as a susceptibility. We recall that in macroscopic electrostatics the susceptibility tensor x relates 
the electric polarization density of a Coulomb system, confined to a region A in a vacuum, to 
the applied electric field. The laws of macroscopic electrostatics allow the components of x to 
be expressed in terms of the dielectric constant of the system. In particular, for a conducting 
ellipse, this theory gives [26] 

^ 27r ^ ' 

where is the length and Ly the width of the ellipse. 

On the other hand, the susceptibility can be related to the microscopic quantities by linear 
response theory. With the electric field applied in the x-direction, this approach gives [27] 

X.. = {{pD - {P.?) (4.2) 
where Px is the x-component of the instantaneous polarization (or equivalently, dipole moment) 

P'-T^jWl- 

Comparing (4.1) and (4.2), and noting that |A| = jL^Ly we obtain 

_ (p^)2 = i-L,(L^ + Ly). (4.3) 
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Now we construct an interval of length 2 in the x-direction from the ellipse by setting = 2, 
Ly = 0. This gives 

iP') - {P.? = ^ (4.4) 
which is in precise agreement with the result (2.5) for the variance of the dipole moment statistic. 
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FIGURE CAPTIONS 



Figure 1. Empirical and theoretical p.d.f. for the statistic J2j=i^j: where the Xj are scaled 

eigenvalues of 2 x 2 GOE matrices. 

Figure 2. Empirical and theoretical p.d.f. for the statistic J2f=i (^^j ~ (^)) where the Xj are 
scaled eigenvalues for the matrix X^X. Here X is a 3 x 2 real rectangular random matrix with 
Gaussian entries. 

Figure 3. Empirical and theoretical p.d.f. for the statistic ^ J2f=i [Oj — tt^ , where the 9j 
(0 < 9j < 2tt) are the phases for the eigenvalues of 3 x 3 CUE matrices. 
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